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We derive the distribution of the synchrotron spectral index across NGC 6946 and investigate the 
correlation between the radio continuum (synchrotron) and far-infrared (FIR) emission using the 
KINGFISH Herschel PACS and SPIRE data. The radio-FIR correlation is studied as a function of star 
formation rate, magnetic field strength, radiation field strength, and the total gas surface brightness. 
The synchrotron emission follows both star-forming regions and the so-called magnetic arms present 
in the inter-arm regions. The synchrotron spectral index is steepest along the magnetic arms (a„ 1) , 
while it is fiat in places of giant Hll regions and in the center of the galaxy {a„ ~ 0.6 — 0.7). The map of 
a-n provides an observational evidence for aging and energy loss of cosmic ray electrons propagating in 
the disk of the galaxy. Variations in the synchrotron-FIR correlation across the galaxy are shown to be 
a function of both star formation and magnetic fields. We find that the synchrotron emission correlates 
better with cold rather than with warm dust emission, when the interstellar radiation field is the main 
heating source of dust. The synchrotron-FIR correlation suggests a coupling between the magnetic 
field and the gas density. NGC 6946 shows a power-law behavior between the total (turbulent) magnetic 
field strength B and the star formation rate surface density Esfr with an index of 0.14 (0.16)±0.01. 
This indicates an efficient production of the turbulent magnetic field with the increasing gas turbulence 
expected in actively star forming regions. Moreover, it is suggested that the B-Espr power law index is 
similar for the turbulent and the total fields in normal galaxies, while it is steeper for the turbulent than 
for the total fields in galaxies interacting with the cluster environment. The scale-by-scale analysis of 
the synchrotron-FIR correlation indicates that the ISM affects the propagation of old/diffused cosmic 
ray electrons, resulting in a diffusion coefficient of Do — 4.6 x 10^* cm'^ s~^ for 2.2 GeV CREs. 

Key words, galaxies: individual: NGC 6946 - radio continuum: galaxies - galaxies: magnetic field - 
galaxies: ISM 



ABSTRACT 



1. Introduction 
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The correlation between the radio and far-infrared 
(FIR) emission of galaxies has been shown to be 
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largely invariant over more than 4 orders of mag- 
nitude in luminosity (e.g. lYun et aL 2001ij) ari d 
out to a redshift of z~3 (e.g. Sargent et al.ll2010l ). 
This correlation is conventionally explained by the 
idea that the FIR and radio emission are both being 
driven by the energy input from massive stars, and 
thus star formation. However, this connection is 
complicated by the observation that the FIR emis- 
sion consists of at least two components; one heated 
directly by massive stars (i.e. the 'warm' dust com- 
ponent), and one heated by the diffuse interstellar 
radiation field or ISRF (i.e. the 'c old' dust compo- 
nent) (see e.g. iDraine et al.l 120071 ). which includes 
the emission from the old stellar population (e.g. 
1X^119901: iBendo et al.ll20T(ih . 

Similarly, the radio emission consists of two 
main components; the thermal, free-free emission 
and nonthermal synchrotron emission. A direct con- 
nection between the free-free emission (of ther- 
mal electrons in Hll r egions) and young massive 
stars is expected (e.g. lOsterbrock fc StockhausenI 
Il960 f) . Conversely, the connection between the syn- 
chrotron emission and massive stars is complicated 
by the convection and diffusion of cosmic ray elec- 
trons (CREs) from their place of birth (supernova 
remnants, SNRs) and by the magnetic fields that 
regulate the synchrotron emission in the interstel- 
lar medium (ISM). 

Hence, warm dust emission /thermal radio emis- 
sion can be directly associated with young stars and 
a correlation between warm dust and thermal radio 
emission is not surprising. On the other hand, the 
connection between cold dust emission /nonthermal 
synchrotron emission and massive stars (and thus 
star formation) is less clear. A better correlation of 
the FIR with the thermal than the nonthermal ra- 
dio emission has already been shown in the LMC, 
M 31 and M 33 (iHughes et al.ll2006t iHoernes et~all 



Table 1. General parameters adopted for 
NGC 6946. 



Il998t iTabatabaei et al.ll2007aD . The CREs experi- 
ence various energy losses while interacting with 
matter and magnetic fields in the ISM, causing the 
power law index of their energy distribution to vary. 
Significant variation of the nonthermal spectral in- 
dex was found in M 33 with a flatter synchrotron 
spectrum in regions of massive SF than in the inter- 
arm re gions and the outer disk (jTabatabaei et all 
l2007bD . 

The critical dependence of the synchrotron 
emission on both magnetic fields and CRE prop- 
agation could cause the nonlinearity in the 
synchrotron- FIR cor relation seen globally for 
galaxy samples (e.g. [Nikiai [iMzS) ■ Propagation 
of the CREs can also cause dissimilarities of the 
synchrotron and FIR morphologies particularly on 
small scales. For instance. Galactic SNRs do not 
seem to be well-c orrelated with the FIR emission 
(e.g. ICohen et~al . 1984). Moreover, within a few 
100 pc of th e star-forming Orion n ebula, no correla- 
tion exists (jBoulanger fc Peraultl IT98S) . In nearby 
galaxies, a lack of correlation on small scales has 
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been shown via detailed multi-scale analysis using 
wavelet transforma ti on (e.g. Hughes et al. 2006; 
iDumas et all [20TTh . iMurphv et al.l ([20061, l2008fl 
showed that recent massive SF could reduce these 
dissimilarities due to generation of a new episode of 
CREs, assuming that the FIR emission is attributed 
to dust heating by the same stars. 

Assuming the massive SF is a s the source of 
both F IR a nd synchrotron emissi on. iHelou fc Bicavl 
(|1993f) and lNiklas fc Beckl (|1997t ) considered a cou- 
pling between magnetic field strength and gas den- 
sity as the reason for the tight radio-FIR cor- 
relation in spite of the sensitive dependence of 
the synchrotron emission on the magnetic field. 
A modified versio n of th is model was suggested 
by IHoernes et al.l (|1998f) to explain the correla- 
tion between the cold dust heated by the ISRF 
and the synchrotron emission, whose energy sources 
are independent. The scale at which the corre- 
lation breaks down provides an important con- 
straint on these models, explaining the scale where 
static pressure equilibrium between the gas and 
CREs/magnetic fields holds. 

Within nearby galaxies, variations in the radio- 
FIR correlation have been sho wn to exist by sev- 
eral a u thors fe.g.lGordon et a l. 2004; Murphv et al.l 



20061: iHuehes et al.l l2006t iMurphv et all 12001 



Dumas et al 11201 ih through a change in the q ratio 



(|Helou et al.lll985l see Sect. 7.3 for definition) or in 
the fitted slope. Furthermore, the smallest scale at 
which the radio-FIR correlation hol ds is not the 
same from one galaxy to another ([Hughes et al.l 
l2006HTabatabaei et alJl2007aHDumas et al.ll201li r 



As the variations in the radio-FIR correlation 
are possibly due to a range of different conditions 
such as the star formation rate (SFR), magnetic 
fields, CRE propagation, radiation field and heating 
sources of dust, this correlation can be used as a 
tool to study the unknown interplay between the 
ISM components and SF. These are addressed in 
this paper through a detailed study of the radio- 
FIR correlation in NGC 6946. 

NGC 6946 is one of the largest spiral galax- 

ies on the sky at a distance of 6.8 Mpc 

([Karachentsev et al.ll2000t ). Its low inclination (38*^) 
makes it ideal for mapping various astrophysical 
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properties across the galaxy (Table 1). NGC6946 
shows a multiple spiral structure with an excep- 
tionally bright arm in the north-east of the galaxy. 
Having several bright giant Hii regions, this Sc 
(SABc) galaxy has active SF as well as strong 
magnetic fiel ds as traced by linearly polarized 
observations (jBeck fc Hoerned I1996D . The global 
star-formation rate is ~7.1M0yr~^ (as listed by 
iKennicutt et al.ll2011l assuming a Kroupa IMF 
and a mass range of 0.1-100 Mq). The dynam- 
i cal mass of this galaxy is ~ 1.9 x 10^^ M0 
(jCrosthwaite fc Turneij|2007[ ). This galaxy harbors 
a mild starburst nucleus (e.g. iBall et al. 1985) and 
there is no stron g evidence for AGN activity (e.g. 
iTsai et al.ll200fil). 

iFrick et "all ( 200l[ ) presented the wavelet analy- 
sis of the radio and mid-infrared (ISOCAM LW3) 
emission in NGC 6946. Here we study this corre- 
lation for dust emission at FIR wavelengths with 
Herschel-PACS and SPIRE from the KINGFISH 
project (Key Insights on Nearby Galaxies: a Far- 
Infra red Survey with Herschel, IKennicutt et al.l 
120111" ) and using various approaches. 

The paper is organized as follows. The relevant 
data sets are described in Sect. 2. After deriving 
the maps of the free-free and synchrotron emis- 
sions (Sect. 3), we derive the distribution of the 
synchrotron spectral index in Sect. 4. We map the 
magnetic field strength in Sect. 5. The radio-FIR 
correlation is calculated using various approaches 
i.e. the q-method, classical pixel-by-pixel correla- 
tion, and as a function of spatial scale in Sect. 6. 
We further discuss the correlations versus magnetic 
fields, SFR, radiation field and gas density (Sect. 7). 
Finally, we summarize the results in Sect. 8. 



lAniano et al.l (|201ll) and resampled to a common 
pixel size of 6" (~170pc). 

The radio continuum (RG ) data at 3.5 and 
20 cm are presented in iBecld (jl99l|) and iBecld 
()2007t ). At 3.5 cm, NGG6946 was observed with 
the 100-m Effelsberg telescope of the MPIfRQ. The 
20 cm data were obtained from observations with 
the Very Large Array (VLAQ) corrected for miss- 
ing short spacings using the Effelsberg data at 
20 cm. To trace the ordered magnetic field, the lin- 
early polarized intensit y data at 6 cm presented in 
iBeck fc HoernesI (|1996[ ) were used. The average de- 
gree of polarization is ~30% for the entire galaxy. 

To investigate the connection between the neu- 
tral gas and the magnetic field, we used the to- 
tal gas (HI -|- H2) mass surface density map which 
was derived using the VLA data of the HI- 21 cm 
line (obtained as part of THINGS, Wal ter et al.l 
120081 ) and the IRAM 30-m C0(2-l) data from 
the H ERA CLES survey as de tailed in iBigiel et al.l 
(|2008[ ) and lLerov et"al1 (|2009t ). 

We used the Ha map of NGC 6946 observed 
with the KPNO 0.9 m in a 23' x 23' field of view and 
with 0.69" pixel"^ (resolution of 1.5"), subtracted 
for the continuum ([Ferguson et al.l[T998l ) and fore- 
ground stars. The contribution fr om the [Nil] line 
emissi on was subtracted following IKennicutt et al] 
|20Q8|) . The Ha emission is corrected for attenu- 
ation by the Galactic cirrus using the extinction 
value given bv lSchlegel etTall (|1998h . The Ha map 



has a calibration uncertainty of ~20%. 

The radio and Ha maps were smoothed to 18" 
resolution (using a Gaussian kernel). All the maps 
were normalized to the same grid, geometry, and 
size before comparison. 



2. Data 

Table 2 summarizes the data used in this work. 
NGC 6946 was observed with the Herschel Space 
Observatory as p art o f the KINGFISH project 
([Kennicutt et al.l 20111) and was described in 

lAniano et al"! 
instrument 



detail in iDale et al.l (|2012[ ) and 
2OI2D . Observations with the PAGS 



(jPoglitsch et al.ll2010D were carried out at 70, 100, 



160 /im in the Scan-Map mode. The PAGS images 
were reduced by the Scanamorphos data reduction 
pipeline (Roussel et al. 2010; Roussel 2012), ver- 
sion 12.5. Scanamorphos version 12.5 includes the 
latest PACS calibration available and aims to pre- 
serve the low surface brightness diffuse emission. 
The 250/im map was observed with the SPIRE in- 
strument ( Griffin et al.. ^010) and reduced using 
the HIPE version spire- 5. 0.1894. T he data were 
subtra cted for the sky as detailed in lAniano et al.l 
(l20T2h . 

The 70, 100, 160 and 250/im images were 
convolved from their native PSFs to a Gaussian 
PSF with 18" FWHM using the kernels from 



3. Thermal/nonthermal separation 

Constraints on the origin and propagation of cos- 
mic rays can be achieved by studying the varia- 
tion in the spectral index of the synchrotron emis- 
sion across external galaxies. To determine the 
variation in the nonthermal radio spectral index, 
we separate the thermal and nonthermal compo- 
nents using a thermal radio tracer (TRT) approach 
in which one of the hydrogen recombination lines 
is used a s a template for the free-free emission 
(see e .g. iDickinson et al.l 120031 : iTabatabaei et al.l 
l2007b[ ). For NGC 6946, we use the Ha line emis- 
sion which is the brightest recombination line data 
available. Both the free-free and the Ha emission 
are linearly proportional to the number of ioniz- 
ing photons produced by massive stars, assuming 



^ The 100-m telescope at Effelsberg is operated by 
the Max-Planck-Institut fiir Radioastronomie (MPIfR) 
on behalf of the Max-Planck-Gesellschaft. 

^ The VLA is a facility of the National Radio 
Astronomy Observatory. The NRAO is operated by 
Associated Universities, Inc., under contract with the 
National Science Foundation. 
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that the Ha emitting me dium is o ptically t hick to 
ioniz ing Lyman photons (jOsterbro ck 1989; iRubinI 
Il968l see also Sect. 3.2). However, the observed 
Ha emission may suffer from extinction by dust 
which will lead to an underesti mate of the free- 
free em issio n. Hence, followin g 'Tabataba ei et al] 
(|2007b[) and iTabatabaei fc Be rkhuiiscn (2013), we 
first investigate the dust content of NGC 6946 in 
an attempt to de-redden the observed Ha emis- 
sion. Then we compare our de-reddening method 
with the one based on a combination of the total 
infrar ed (TIR) and Ha luminosity (Kennicutt et al.l 
I2009D . 

3.1. De-reddening of the Ha emission 

Following iDraine fc Lil (|2007l) . the intersteUar 
dust heating has been modeled in NGC 6946 
(jAniano et al.l I2012D assuming a ^-function in ra- 
diation field intensity, U, coupled with a power-law 
distribution Umin < U < J/max, 



dMdust/dU = Mdust (1 



7 



7)SiU- 
a- 1 



Umm) 
■U' 



(1) 



where U is normalized to the local Galactic inter- 
stellar radiation field, Mdust is the total dust mass, 
and (1 — 7) is the portion of the dust heated by 
the diffuse interstellar radiation field defined by 
U = ?7min- The minimum and maximum interstel- 
lar radiation field intensities span 0.01 < t/mi n < 30 
and 3 < logC/^ax < 8 (see iDale et al.l I2OOIL 
l2012f) . Fitting this model to the dust SED covermg 
the wavelength range between 3.5 /im and 250 /im, 
pixel-by-pixel, results in the 18" maps (with 6" pix- 
els) of the dust mass surface density (a), the dis- 
tribution of the radiation fields (U), and the total 
infrared (TIR) luminosity emitted by the dust. 



The map of a is equivalent to a map of the 
dust optical depth at Ha wavelengths given by 
THa = f i^Ha wherc khq is the dust opacity. Taking 
into account both absorption and scattering, Kua = 
2.187 X lO'^cm^g^^, assuming a Milky- Way value 
of the total/selectiv e extinction ratio of Rv=3.1 
(jWeingartner fc Dra inc 2001). The distribution of 
thq over the disk of NGC 6946 overlaid with con- 
tours of the Ha emission is shown in Fig. 1. We 
note that the dust mass from modeling the SED 
at 18" resolution may be overestimated by ~20% 
due to the lack of longer wavelength constraints 
(jAniano et al.ll2012l) . 

Regions with considerable dust opacity (thq > 
0.7) follow narrow dust lanes along the spiral arms 
(e.g. the narrow lane in the inner edge of the bright 
optical arm) and are found mainly in the inner 
disk. High opacity dust is found in the center of 
the galaxy with thq — 5. This corresponds to 
a silicate optica l depth of T9 . 7 ::^0 .5 which is in 
agreement with ISmith et al.l (j2007t) . In the cen- 
tral '^60pc which is much smaller than our resolu- 
tion, much larger estimates of th e extinction have 
been found using total gas masses (jSchinnerer et al.) 
MM)- 

Figure 1 also shows that the Hii complexes are 
dustier in the inner disk (with thq ^1-4) than 
in the outer parts (with trq <0.6) of NGC 6946. 
Across the galaxy, the mean value of thq is 0.43 ± 
0.04 (median of 0.34±0.04). Therefore, apart from 
the center, NGC 6946 is almost transparent to pho- 
tons with A ~ 6563A propagating towards us. 

The thq derived can then be used to correct the 
Ha emission for attenuation by dust, taking into ac- 
count the effective fraction of dust actually absorb- 
ing the Ha photons. Since these photons are usually 
emitted from sources within the galaxy, the total 
dust thic kness Tun only provide s an upper limit. 
Following iDickinson et al.l (j2003l) , we set the effec- 
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RIGHTASCENSION tJZJOO) RIGHTASCEN3IQN (J2000) 



Fig. 1. Left: Dust optical depth at Ha wavelength tho overlaid with contours of the Ha emission from 
NGC 6946. Contour levels are 1, 3, 6, 15, 30, 40 ^erg s~^ cm^^ sr^^. The bar at the right shows the values 
of THa ■ Right: Radiation field U map (mass-weig hted mean starlight heating intensity) of NGC 6946 in 
units of radiation field in solar neighborhood Uq (jAniano et al.|[201^ . indicated by the bar at the right 
of the image. The beam area is shown in the lower left corners. 



tive thickness to Toff — fd ^ thq with fd being the 
dust fraction actually absorbing the Ha; the atten- 
uation factor for the Ha flux is then e~'^''". At our 
resolution of 18" ~530pc, one may assume that the 
Ha emitting ionized gas is uniformly mixed with the 
dust, which would imply fd — 0.5. Considering the 
fact that the ionized gas has a larger extent than 
the dust, Dickinson et al. (2003) found a smaller 
effective factor {fd = 0.33) based on a cosecant law 
modeling. We also adopt fd = 0.33 for NGC 6946. 
We note that this choice barely influences the ther- 
mal fraction of the radio emission, due to the small 
THa (Sect. 3.2). 

Of course, it would be preferable not to use 
a uniform value fd for the whole galaxy, but one 
that is adapted to the geometry (w ell mixed diffuse 
medi um or shell-like in Hii regions, iWitt fc Gordm] 
l2000f ) and the dust column density. However, this 
would require specifying the location of the stel- 
lar sources and the absorbing dust along the line 
of sight and solving the radiative transfer problem 
with massive numerical computations which is far 
from our first order approximation. 

In another approach, assuming that the dust 
is mainly heated by the massive stars, we cor- 
rected the Ha emission by combining it with the 
TIR (integrated dust luminosity in the 8-1000 /im 
wavelength range): Hacoir= Haobs + 0.0024 TIR 



(jKennicutt et al.l2009l) . Interestingly, this approach 
is linearly correlated with the de-reddening using 
Toff (Fig. 2), with an offset of 0.19 dex and a disper- 
sion of 0.14 dex. Figure 2 shows that at the highest 
luminosities, the corrected Ha values agree, corre- 
sponding to the calibration of the second approach 
specifically to star-forming regions. Outside of these 
regions, the Ha-TIR ratio approach overestimates 
the correction applied to the observed Ha, probably 
because of contributions from other dust-heating 
sources. Masking out the diffuse emission in the 
inter-arm regions and outer disk (i.e., considering 
only the spiral arms and SF regions) , both the offset 
and dispersion reduces to 0.11 dex and hence both 
methods agree within the uncertainties (~ 20% due 
to calibration). This likely indicates that the diffuse 
dust is not heated by the UV radiation of ionizing 
stars. 



3.2. Distribution ofttie free-free and syncftrotron 
emission 

Using the corrected Ha emission from the first ap- 
proach, we derive the intrinsic Ha intensity, Iq, 
according to / ~ /q e~'^'=^^ . Integration of the de- 
reddened Ha map out to a radius of 11.9 kpc (414") 
yields a luminosity of Lhq — (3.46 ± 0.05) x 
lO^^ergs"^ that is higher than the foreground- 
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Fig. 3. Free-free (lejf) and synchrotron emission (right) from NGC6946 at 3.5 cm (top) and 20 cm {bot- 
tom). The angular resolution is 18" (shown in the lower left corner of the panels) with a grid size of 6". 
The bars at the right of the images show the intensity values in mJy/beam. Note that the areas are not 
the same at 20 cm and 3.5 cm due to their different ob served fields. At 20 cm the synchrotron emission is 
overlaid with contours of linearly polarized intensity (|Beck fc HoernedllQQGl ). The contour levels are 70, 
120, 160/uJy/beam. 



corrected luminosity by ~20%. A similar increase 
in the Ha flux has been derived in other nearby 
galaxies, M33 (~ 13%, 'Tabatabaei ct al. 200'^ 
and M31 (~30%, .Tabatabaei fc BerkhuiisenI 

120101). 

[Dickinson et al.l ()2003[ ) showed that the Ha 
emitting medium in our Galaxy is optically thick 
to io nizing Lyman photons (case B, lOsterbrockl 
119890 not only for Hii regions (rLya ~ 10^ - 10^°) 
but also for faint Ha features at intermediate and 
high Galactic latitudes (rLya 1 — 30). Assuming 
the same condition applies for NGC 6946, the emis- 
sion measure (EM) follows from the Ha intensity 

in units of ergcm~ ^s~^ sr~^) via the expression 

Valls-Gabaudlll998h : 

/Ha = 9.41 X 10"^Tg4^ °^^10~^ EM , (2) 



where the electron temperature, Te4, is in units of 
10'' K, and EM in cm~^ pc. The emission measure 
is related to the continuum optical thickness, Tc, of 
the ionized gas by 

= 3.278 X 10"'^aT^4^-^^i^GHz (1 + 0.08) EM , (3) 

with a ~ 1 (|Dickinson et al.l I2003D . The factor 
(1 -I- 0.08) takes into account the contribution from 
singly ionized He. The brightness temperature of 
the radio continuum (free- free) emission, T{,, then 
follows from 

n^Te{l~e-^^) . (4) 
Eq. (4) with Eqs.(2) and (3) gives: 

rTfc = Te(l-e~^^H.)^ ^^^^ 
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Fig. 2. Ha luminosity de- reddened based on Tcs 
versus de-reddening using the TIR luminosity. Also 
shown are the lines of 1:1 correspondence (dashed) 
and Y=X - 0.19 (solid). 



Hence, the free-free emission can be derived 
separately at each radio wavelength. The result- 
ing distributions of the intensity of the free-free 
emission in mJy/beam at 3.5 and 20 cm are shown 
in (Fig. 3, left panels) foiQ = 10^ K. Using a 
constant electron temperature is supported by the 
shallow metallicity gra dient found in this galaxy 
(jMoustakas et al.l l2010D . Subtracting the free- free 
emission from the observed radio continuum emis- 
sion results in a map of the synchrotron emission 
(Fig. 3, right panels). 

The synchrotron maps exhibit diffuse emission 
extending to large radii indicating diffusion and 
propagation of the CREs. Strong synchrotron emis- 
sion emerges from the galaxy center, giant star- 
forming regions, and spiral arms, which could be 
due to stronger magnetic fields and/or young- 
energetic CREs close to the star-forming regions. 
Interestingly, the so called 'magnetic arms' traced 
by th e linearly polarized intensity (Beck fc Hoerne^ 
Il996f ) are clearly visible in the 20 cm synchrotron 
map. The fact that they are less prominent at 
3.5 cm implies that these arms are filled by older 
and lower energetic CREs (see Sect 4.2). The ther- 
mal free-free map, on the other hand, exhibits nar- 
row spiral arms dominated by the star-forming re- 
gions. 

Integrating the observed, synchrotron and free- 
free maps in the plane of the galaxy (i=38°) around 
the center out to a radius of 324"(9.2 kpc), we ob- 
tain the total flux densities and thermal fractions at 
3.5 and 20cm (Table 3). The thermal fractions are 
about 18% and 7% at 3.5 and 20 cm, respectively. 
As mentioned before, we assumed a dust attenua- 



^ Smaller values of Te are repo rted from the mea- 
surements in the M ilky Way (e.g. iHaffner et al.|[l999l : 
iMadsen et aD 120061 '). Assuming Te = 7000 K, the ther- 
mal fraction would decrease by about 23%. 



Table 3. Global radio continuum fiux densities and 
thermal fractions in NGC 6946. 



A 


Observed 


Free- free 


Thermal 


(cm) 


flux density 
(mJy) 


flux density 
(mJy) 


fraction 

% 


3.5 


422 ± 65 


78 ± 10 


18.4 ±3.7 


20 


1444 ± 215 


97 ± 13 


6.7±1.3 



tion factor of fd = 0.3. For a uniform distribution 
of dust and ionized gas (fd — 0.5), the thermal 
fractions increase to about 21% and 8% at 3.5 and 
20 cm, respectively. 

4. Synchrotron spectral index 

Using the nonthermal radio fiuxes at 3.5 and 20 cm, 
we obtained the spectral index of the nonthermal 
radio emission. This was only computed for pixels 
with flux densities of at least three times the rms 
noise at both frequencies. The synchrotron spectral 
index, a„, shows a smooth variation across the disk 
of NGC 6946 (Fig. 4). The greenish color in Fig. 4 
corresponds to the synchrotron emission with flat 
spectrum (a„ < 0.7), and the reddish color to the 
regions with a steep spectrum (q;„ > 0.9), i.e. emis- 
sion from lower-energy CREs. The flat spectrum is 
found in giant Hll regions and the steep spectrum 
in the inter-arm regions, the eastern part of the cen- 
tral galaxy, and the south of the major axis where 
an = l.OibO.l. As shown in the histogram represen- 
tation (Fig. 5), the median value of a„ across the 
galaxy is 0.81 with a dispersion of 0.18. Figure 5 
also presents a histogram of the 'observed' spectral 
index a (obtained using the observed radio data at 
3.5 and 20cm, i.e., contaminated by the free-free 
emission), for comparison, with a ~ 0.7 on the av- 
erage and with a dispersion of 0.16. 

4.1. Synchrotron spectral index versus star formation 

In the star-forming regions, the synchrotron spec- 
trum is relatively flat with an average index of 
a„=0.65±0.10, the typical spectral index of young 
CREs in SF regions (in supernova remnants, a„ is 
about 0.5 on the avera ge and could even be flatter, 
see e.g. [Gordon et al.l [19991 iRevnolds etlll 120121 : 
lLongaiilll994l and references therein). Table 4 lists 
an in 9 giant Hii regions annotated in Fig. 4, (called 
Enuc follo wing the nomenclatu re and location pre- 
sented in iMurphv et al.ll20To[ ). Enuc2 and Enuc3 
have a relatively steep spectrum which could be 
due to an energy loss of CREs in a magneto-ionic 
medium along the line of sight seen at the position 
of these Hii regions. This is possible for Enuc3 since 
this source is adjacent to the strong northern mag- 
netic arm (see below). In the neighborhood of the 
Enuc2, however, no strong polarized emission is de- 
tected. The estimated spectral index of this source 
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Table 4. Thermal fractions and synchrotron spec- 
tral index of the giant Hii regions shown in Fig. 4. 



Object 


7th 

(%) 


/th 

(%) 


an 


Enucl 


50 ±2 


26 ±1 


0.66 ±0.03 


Enuc2 


80±4 


41 ±2 


0.90 ±0.08 


Enuc3 


73 ±3 


34 ±2 


0.86 ±0.07 


Enuc4 


41 ±2 


26 ±1 


0.48 ±0.05 


Enuc5 


68 ±1 


34 ±1 


0.75 ±0.05 


Enuc6 


49 ±3 


23 ±2 


0.70 ±0.03 


Enuc7 


56 ±2 


27 ±2 


0.74 ±0.04 


Enuc8 


47 ±2 


20 ±1 


0.54 ±0.02 


Enuc9 


48 ±5 


23 ±3 


0.70 ±0.03 



Fig. 4. Synchrotron spectral index map of 
NGC 6946 overlaid with contours of the lin- 
early polarized emission at 6 cm (top) tracing the 
ordered magnetic field in the sky plane with the 
same levels as in Fig. 3 and of bright Ha sources 
{bottom). The Ha contour levels are 18, 24, 35, 
and 47% of the maximum intensity. The 9 giant 
Hii regions are indicated (Table 4). The bars at 
the right of the images show the values of the 
synchrotron spectral index. 



could be affected by an underestimate of the ob- 
served 3.5 cm emission, since this source sits on the 
edge of the 3.5 cm observed field (primary beam). 

Considering synchrotron, inverse-Compton, ion- 
ization, and bremsstrahlung cooling mechanisms of 
CREs as well as usin g a prescription for t he es- 
cape of these particles, iMurphv et al.l (|201Cl[ ) mod- 
eled the radio SEDs and found an a^ — 0.8 for 



the giant Hii regions in Table 4, which is steeper 
than our finding and equals to the average a„ 
for the entire galaxy. One contributing factor to 
their steeper indices could be a result of assum- 
ing a fixed ISM density for the extranuclear re- 
gions of 0.1 cm^'^ (motivated by their beam size 
area of 30" ^ 0.9 kpc), leading to dominant syn- 
chrotron and IC losses rather than the ionization 
and bremsstrahlung cooling mechanisms of CREs 
(and hence a rather steep spectrum for these ob- 
jects). In the present work, we determine a„ at 
a smaller beam size of ^0.5 k pc (the scale of the 
giant SF regions in NGC 6946 iKennicutt fc EvansI 
I2012D and without any assumption on the ISM den- 
sity and/or cooling mechanism of CREs. On these 
scales, more energetic CREs and/or stronger mag- 
netic fields close to the SF regions provide the syn- 
chrotron emission with a flatter a„ on the average. 

Table 4 also lists the thermal fractions at 3.5 cm 
and 20 cm in these extranuclear Hli complexes. The 
thermal fractions at 3.5 cm are higher than those at 
20 cm by a factor of about two (apart from Enuc4) . 
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4.2. Synchrotron spectral index versus magnetic 
fields 

The magnetic fields i n NGC 6946 hav e bee n 
extensively s tudied by iBeck fc HoernesI (|1996[ ). 
iRohde et all (|1999[ ). and iBec^ (|2007D. The con- 
tours in Fig. 4 show the linearly polarized inten- 
sity (PI) at 6.3 cm which determines the strength 
of the order ed magnetic field in th e plane of the 
sky (see e.g. iTabatabaei et al.ll2008h . Interestingly, 
there is a good correspondence between the steep 
synchrotron emission and th e PI, particularly along 
the northern magnetic arm (jBeck fc Hoerneslll996[ ) 
and also along the strong ordered magnetic field 
in the central disk (anisotropic turbulent magnetic 
field, iBeckl 1200711 . In these regions, the spectral 
index of = 1.0 ± 0.1 indicates that CREs 
suffer strong synchrotron losses propagating along 
NGC 6946's ordered magnetic field. 

Overall, the synchrotron spectral index map 
agrees with the energy loss theory of relativistic 
electrons propagating away from their origin in 
star- forming regions in the ISM (e.g. see Chapter 
18 of lLongaiij[T99llBiermann et al.ll200H) . The dif- 
ference in a„ in star-forming regions and along the 
magnetic arms in NGC 6946 is similar to that pre- 
dicted bv lFletcher et al.| ()2011l) for the spiral galaxy 
M51. 



5. Maps of total and turbulent magnetic 
fields 

The strength of the total magnetic field Btot can 
be derived from the total synchrotron intensity. 
Assuming equipartition between the energy densi- 
ties of the magnetic field and cosmic rays (ecr = 

Btot = C(a„,X,i)[/„]^, (6) 

where /„ is the nonthermal intensity and C is a 
function of q;„, K the ratio between the num- 
ber densities of cosmic ray protons and elec- 
trons, and L the pathlength th rough the syn- 
chrotron emitting medium (see iBeck fc Krausd 
120051 : ITabatabaei et al.ll2008^ . Using the maps of 
/„ and a„ obtained from the TRT method and as- 
suming that the magnetic field is parallel to the 
plane of the galaxy (inclination of i — 38° and po- 
sition angle of the major axis of PA=242°), Btot 
is derived across the gala xy. In our calculation s, we 
apply values of ~ 100 (|Beck fc Krausell2005^ and 
L ~ Ikpc/cosz. Figure 6 shows strong Btot in the 
central region of the galaxy, the arms and the star- 
forming regions. 

As a fraction of the polarized intensity PI is re- 
lated to the strength of the ordered magnetic field, 
and the nonthermal intensity /„ to the total mag- 
netic field in the plane of the sky, /„ — (PI/0.75) 
gives the nonthermal emission due to the turbu- 
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Fig. 6. Strength of the total Btot (top) and tur- 
bulent Btur (bottom) magnetic fields in NGC 6946. 
The bars at the top of each image show the mag- 
netic field strength in /i Gauss. 

lent magnetic field Btufl Using this intensity with 
Eq. dS]) yields the distribution of Btur across the 
galaxy. Similar to Btot, Btur is higher in places of 
star-forming regions (Fig. 6, bottom). We note that 
any ordered field which is not resolved and depo- 
larized within the beam would contribute to the 
turbulent magnetic field. In other words, we cannot 
distinguish between unresolved structures of the or- 
dered field and truly turbulent fields for structures 
smaller than the beam. 

The magnetic field strength estimated using the 
'standard method' for thermal/nonthermal separa- 



* In a completely ordered magnetic field, the maxi- 
mum degree of linear polarization is about 0.75 (e.g. 
iWestfold,1959i ) 
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tion, i.e. assuming a fixed synchrotron spectral in- 
dex (see Appendix A), is similar to the TRT es- 
timate (within 10% difference) in the spiral arms, 
and is underestimated by 35% in the giant Hii re- 
gions, and 45% in the nucleus. This is because, in 
those regions, /„ is smaller and q;„ is larger based 
on the standard method. 



6. Radio-FIR Correlation 

The radio correlations with both monochro- 
matic and bolometric FIR observations are ob- 
tained using three different approaches: classical 
pixel-t o-pixel (Pearson c orrelation), wavelet multi- 
scale dFr ick et alj l200lD . and the q-ratio method 
(|Helou e t al. 1983) ■ Although the main goal of this 
paper is to investigate the synchrotron~FIR corre- 
lation, we also present the correlations between the 
emissions of the thermal free-free and each of the 
monochromatic FIR bands. 



6.1. Classical correlation 

The pixel-to-pixel correlation is the simplest mea- 
sure of the correlation between two images, /i (x, y) 
and f2{x,y), with the same angular resolution and 
the same number of pixels using Pearson's linear 
correlation coefficient, rj£|: 



v/S(/l.-(/l))2S(/2*-(/2))2 



(7) 



When the two images are identical, Tj, = 1. The 
correlation coefficient is rc=-l for images that are 
perfectly anti-correlated. The formal error on the 
correlation coefficient depends on the strength of 
the correlation and the number of independent pix- 
els, n, in an image: Arc = a/I — r^/^n — 2. 

We calculated the correlations between both the 
radio free- free/synchrotron and the FIR 70, 100, 
160, and 250/im maps, restricting the intensities 
to > 3 X rms noise. We obtained sets of indepen- 
dent data points (n) i.e. a beam area overlap of 
< 20%, by choosing pixels spaced by more than 
one beamwidth. Since the correlated variables do 
not directly depend on each other , we fitted a power 
law to the bisector in each case (jlsobe et al.lll990l: 
iHoernes et al.l[T998h . The Student's t-test is also 
calculated to indicate the statistical significance 
of the fit. For a number of independent points of 
n > 100, the fit is s ignificant at the 3a level if 
t > 3 (e.g. IWallll 19791) . Errors in the slope b of the 
bisector are standard deviations (Ict). 

The results for both radio wavelengths are pre- 
sented in Table 5. The calculated Student's t-test 

^ Please note that using this method it is not possi- 
ble to separate variations due to a change of the physi- 
cal properties with scale e.g. propagation of CREs (see 
Sect. 6.2). 



values are large (> 33) indicating that the fitted 
slopes are statistically significant. With coefficients 
of Tc > 0.8, good correlations hold between the 
FIR bands and observed radio (RC)/free-free emis- 
sion. The FIR correlation coefficients with the syn- 
chrotron emission are slightly lower than those with 
the free-free emission. 

The synchrotron emission is slightly better cor- 
related with the 250^m emission (as a proxy for 
cold dust) than with the 70/im emission (as a 
proxy for warm dust). On the contrary, a free- free- 
cold/ warm dust differentiation is only hinted, but 
not yet clear given the errors from the Tc values. 

The free-free emission exhibits an almost linear 
correlation with the warmer dust emission at 70 
and 100/im (with a slope of 6 ~ 0.9). The correla- 
tion becomes more and more sub-linear with dust 
emission probed at 160 and 250/im. 

The synchrotron emission, on the other hand, 
tends to show a linear correlation with colder rather 
than warmer dust (with a super-linear correlation) . 
A similar trend is also seen between the observed 
RC and FIR bands. 

Super-linear radio-FIR correlati ons have been 
also fo und f or samples of g alaxies bv I Price fc Durid 
(|19921) and iNiklasI (|l997ali . which were attributed 
to the non-linearity of the synchrotron-FIR corre- 
lation and/or to the fact that colder dust may not 
be necessarily heated by the young massive stars. A 
better synchrotron-cold than -warm dust correla - 
tion was also found in M31 (jHoernes et a l.' iggS*) 
and in a sample of late- type galaxies ( Xu et al.l 
IT99I . 

One possible issue is that our use of the dust 
mass to de-redden the Ha emission, and thus the 
free-free emission, is somehow influencing the corre- 
lations. To test this, we re-derive the free-free emis- 
sion using the observed Ha emission (not corrected 
for extinction) and re-visit the thermal/nonthermal 
correlations with the FIR bands. The results are 
given in parenthesis in Table 5. The differences are 
less than 4% and within the errors. 

The synchrotron emission based on the stan- 
dard method also shows a decrease of the slope 
with increasing FIR wavelength (Table 6), similar 
to the TRT based study. However, the expected 
better linearity of the free-free-warmer dust is not 
seen using the standard thermal/nonthermal sepa- 
ration method. As shown in Sect. 5, in our resolved 
study, this method results in an excess of the free- 
free diffuse emission in the inter-arm regions where 
there is no warm dust and TIR counterparts. This 
is most probably caused by neglecting variations of 
a„ locally across the galaxy, since in global studies, 
the standard separation method le ads to a linea r 
thermal radio-FIR correlation (e.g. lNiklasl[T997bl) . 
Based on the standard method, the separated RC 
components are not as tightly correlated with the 
FIR bands as the observed RC-FIR correlations. 
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Table 5. Radio-IR relations on pixel-by-pixel basis in N6946. 



X 


Y 


b 


' c 


t 


]2 


RC(3.5) 


1 70 


1.49±0.04 


0.90±0.02 


53.0 


662 




IlOO 


1.49±0.03 


0.92±0.01 


61.8 


696 




I16O 


1.20±0.03 


0.88±0.01 


55.0 


883 




T 

i250 


1 no_i_n no 


n on_Ln no 


06.3 


834 


FF(3.5) 


T 

iro 


U.92(U.yoj±(J.U2 


U.yU(U.88j±U.U2 


52.5 


650 




Iioo 


n RR{() Ql Wo 09 




00.0 






IlGO 


fl 77^0 7q\+c\ m 


88^0 88W0 02 


51.3 


770 




I250 


n fi7(n fiqWn ni 

U.U 1 lU.Ut? J _I_U . U J- 


87f0 8fiWo 09 


48.2 


750 


SYN(3.5) 


1 70 


1.60(1. 58)±0.04 


0.80(0.82)±0.02 


33.8 


644 




Iioo 


1.59(1. 57)±0.04 


0.84(0.84)±0.02 


40.2 


677 




I16O 


1 45^ 43Wn 03 


8fif0 87W0 09 


47.1 


784 




I25O 


1 1 Sn 1 7Wn 03 

-L . -L 1 J- . -L 1 J _I_U . UO 


8fif0 87^+0 09 

U.OUIU.O/ 1 _I_U .u^ 


46.4 


761 


RC(20) 


1 70 


1.52±0.03 


0.87±0.02 


45.4 


664 




Iioo 


1.48±0.02 


0.90±0.02 


54.4 


697 




I16O 


1.16±0.02 


0.92±0.01 


70.1 


911 




T 

-1-250 


i.UZzbU.Ui 


n no 1 n n 1 

u.y2±u.ui 


68.3 


850 


FF(20) 


1 70 


0.92dz0.02 


0.90±0.02 


50.3 


656 




Iioo 


U.o ( rnU.UZ 


QO-i-O 09 


f^9 ^ 


Ooo 




I16O 


0.79±0.01 


0.88±0.02 


52.3 


800 




I25O 


0.68±0.01 


0.87±0.02 


49.3 


793 


SYN(20) 


1 70 


1.53±0.03 


0.83±0.02 


38.2 


661 




Iioo 


1.48±0.03 


0.87±0.02 


65.8 


695 




I16O 


1.19±0.02 


0.91±0.01 


65.7 


900 




I25O 


1.04±0.02 


0.90±0.01 


60.0 


846 



Notes. Pearson correlation coefficients Vc and bisector slopes b, log(Y)oc 61og(X), between the total RC, Free- 
Free (FF), synchrotron (SYN) from the TRT method and the far-infrared PACS and SPIRE bands. The values in 
parenthesis show the corresponding correlations when the observed Ha emission is directly used as the free-free 
tracer. Ordinary least-squares fits of bisector log(Y)~ blog{X) are given through n pairs of (logX, logY), where 
n is the number of independent points (Isobe et al. 1990); t is the Student's t-test. 



Table 6. Same as Table 5 for the 3.5 cm Free-Free 
(FF)/synchrotron (SYN) emission separated using 
the standard method (see Appendix A.). 



X 



FF(3.5) 



Y 



1 70 

Iioo 

I16O 
I25O 



SYN(3.5) I70 
Iioo 
I160 
I250 



1.17±0.03 
1.13±0.03 
1.03±0.03 
0.91±0.03 

1.60±0.05 
1.51±0.04 
1.28±0.03 
1.13±0.03 



0.78±0.03 
0.75±0.03 
0.72±0.03 
0.72±0.03 

0.77±0.03 
0.79±0.02 
0.84±0.02 
0.83±0.02 



28.4 
25.5 
24.3 
23.8 

29.7 
33.9 
38.1 
37.5 



508 
518 
543 
540 

594 
605 
626 
622 



emission with a large-scale distribution, typically. 
However, such a correlation can be misleading when 
a bright, extended central region or an extended 
disk exists in the galactic image. In this case, the 
-classical correlation is dominated by the large-scale 
^structure, e.g. the disc of the galaxy, while the 
(more interesting) correlation on the scale of the 
spiral arms can be much worse. The classical cor- 
relation provides little information in the case of 
an an ti-correlation on a specific scale ()Frick et al.l 
l200lh . U sing the 'wavelet c ross-correlation' intro- 
duced bv lFrick et all (|200lD . one can calculate the 
correlation between different emissions as a func- 
tion of the angular scale of the emitting regions. 
The wavelet coefficients of a 2D continuous wavelet 
transform are given by: 



6.2. Multi-scale correlation 

Classical cross-correlations contain all scales that 
exist in a distribution. For example, the high- 
intensity points represent high-emission peaks on 
small spatial scales belonging to bright sources, 
whereas low-intensity points represent weak diffuse 



-)rfx', (8) 



VF(a,x) 



where /(x) is a two-dimensional function (an im- 
age), V'(x) is the analyzing wavelet, the symbol * 
denotes the complex conjugate, x — {x,y) defines 
the position of the wavelet and a defines its scale, k 
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is a normalization parameter (we use k = 2, which 
is equivalen t to e nergy norm alization) . F ol lowin g 
iFrick et alj (|200lD and, e.g., iLaine et all (|2010t ). 
we use the 'Pet Hat ' function as the analyzing 
wavelet ip{x.) to have both a good scale resolution 
and a good spatial resolution. The wavelet cross- 
correlation coefBcient at scale a is defined as 



.(«) 



[Afi(a)Af2(a)]i/2 ' 



(9) 



where the subscripts refer to two images of the 
same size and linear resolution and M{a) = 
\W{a, x)pdx is the wavelet equivalent of 
the power spectrum in Fourier space. The value of 

varies between -1 (total anti-correlation) and 
+1 (total correlation). The correlation coefficient of 
l^'tol = 0.5 is translated as a marginal value for ac- 
ceptance of the correlation between the structures 
at a given scale. The error in is estimated by 
the degree of correlation and the number of inde- 
pendent points (n) as Ar^{a) = ^1 — T\,l\Jn — 2, 
where n ^ LxLy/a?, and {Lx,Ly) are the map 
sizes in {x,y). Thus, towards larger scales, n de- 
creases and the errors increase. 

To prevent a strong influence of the nucleus 
on the wavelet analysis, the central 2 kpc was sub- 
tracted from the PACS and SPIRE maps as well as 
the radio maps of NGC6946. Then, the maps were 
decomposed using Eq. (8) into 12 angular scales (a) 
between 32"(-0.9kpc) and 313"(-9kpcB The de- 
rived wavelet coefficients W{a,x) of the radio and 
FIR maps were then used in Eq. (9) in order to 
derive the radio~FIR correlation coefficient on 
the 12 scales of decomposition. Figure 7 shows r^, 
between the maps of the FIR bands and the 20 cm 
synchrotron, free-free, and observed RC versus spa- 
tial scale a. 

The total RC-FIR and the synchrotron-FIR 
correlations are higher for scales > 4 — 5 kpc, the 
scale of the diffuse central disk. On smaller scales 
the synchrotron-FIR correlations fluctuate (within 
the errors) close to the threshold value r^=0.5. 

On the smallest scale of 0.9 kpc, correspond- 
ing to the width of the complexes of giant molec- 
ular clouds (GMCs) and star formation within 
spiral arms (typical widt h of spiral arms in the 
galaxy- 1.6 - 2 kpc, iFrick et al.l l200lD . the 
synchrotron-colder dust correlation (e.g. 250/im, 
r.u] — 0.53±0.02) is slightly better than the 
synchrotron-warmer dust one (70yum, = 
0.45±0.02). Right after a small maximum on scales 
of spiral arms {'^ 1.6 — 2 kpc), a small minimum 
in the RC/synchrotron-FIR correlation occurs at 
a ~ 2.5 kpc, the scale covering the width of the spi- 
ral arms plus diffuse emission in between the arms. 
Hence, the minimum could be due to the diffuse 



® The wav elet analysis was per formed using the code 
described bv lDumas et al.l (|2011l '). 



nature of the nonthermal emission caused by the 
propagation of CREs. The existence of the 'mag- 
netic arms' where there is no significant FIR emis- 
sion, could also cause the general weaker correlation 
on scales of their widths, i.e., < 3 kpc (see Sect. 7.7). 

The FIR bands are better correlated with 
free-free than with synchrotron emission on a < 
3 kpc, as expected (note that this difference is not 
that obvious in the classical correlation at 20 cm). 
Moreover correlations show a split in terms of the 
FIR band, with warmer dust having a better cor- 
relation than colder dust. Such a split is again not 
obtained from the classical correlation. This occurs 
because the classical correlation is biased towards 
large scales (due to the bright disk of NGC6946), 
where the cold/ warm dust split is less pronounced 
(or disappears in the free-free-FIR plot , Fig. 7 mid- 
dle). 

6.3. FIR/ radio ratio map 

In order to determine variations in the radio-FIR 
correlation across the disk of NGC 6946, we con- 
structed a logarithmic FIR/radio ratio map using 
the 20 cm 'sy nchrotron ' map and a FIR42_i22/jm 
(following Helo u et al.l [198 luminosity map de- 
rived by integrating the sp ectral energy distri- 
bution (SED) pixel-by-pixel (lAniano et aLl 



According to the convention of iHelou et al.l 
q is defined as: 



q = log ( 



FIR 



3.75 X 10i2Wm- 



— loe 



^Hz 



-I 



with 5*1.4 the synchrotron 20 cm flux. The resulting 
map is shown in Fig. 8. Only pixels detected above 
the 3cr level in each radio and FIR map were con- 
sidered. The mean q across the galaxy is 2.22±0.10 
with 0.10 the standard deviation. In star- forming 
regions we find q = 2.30-2.40 whi ch is in good agree - 
ment with the value derived by lYun et al] ()2001b[ ) 
for a sample of nearby galaxies (with a mean value 
oiq — 2.34). In the inter-arm regions, q decreases to 
values lower than 2.0. Hence, spiral arms are easily 
discernible in the q map. 

The free-free emission provides a proper mea- 
sure for the present -day star formation rate (e.g. 
iMurphv et"aLl l201lD . The upper panel in Fig. 8 
shows a good correlation between q and the free- free 
emission (contours), in agreement with the impor- 
tant role of SF in controlling the radio-FIR correla- 
tiorfl. Besides SF, magnetic fields also seem to play 
a role as suggested by Fig. 8 (middle) where we ob- 
serve a correspondence between small q values and 
peaks of the polarized intensity. Furthermore, there 



^ We stress that the arm/inter-arm variation of q as 
well as the dependence of q on SFR is a direct conse- 
quence of the nonlinearity of the radio-FIR correlation. 
Hence, in general, averaging q values could be mislead- 
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is a good agreement between high q values and the 
turbulent magnetic field (Fig. 8 bottom). Is there 
any dependency or competition between SF and 
magnetic fields in controlling the radio-FIR corre- 
lation? A more quantitative comparison between q, 
SFR, and magnetic fields is given in Sect. 7. 

7. Discussion 

In the following, we discuss the variations in 
the synchrotron-FIR correlation across the disk of 
NGC 6946 as a function of radiation field (impor- 
tant for the dust heating), star formation and mag- 
netic fields. We further investigate possible connec- 
tions between SF and the ISM components, mag- 
netic fields and gas. This will help us to understand- 
ing the origin of the observed synchrotron-FIR cor- 
relation in this galaxy. 

7.1. Correlations versus radiation field 

The spatial variations in q may mean that the slope 
of the correlation is not constant. The higher q, 
e.g., in the arms rather than in the inter-arm re- 
gions can be translated as a larger slope (6) of 
the synchrotron-FIR correlation in the arms than 
in the inter-arm regions (a l arger arm/inter-arm 
slope was also found in M51 bv lDumas et al.|[2011[ ). 
Similarly, the good correspondence between q and 
star formation (Sect. 7.3) implies that the slope of 
the synchrotron-FIR correlation is higher in SF re- 
gions and lower in between the arms and outer disk. 
However, if we consider that the synchrotron-FIR 
correlation originates from massive SF, no corre- 
lation would be expected in regions devoid of SF. 
Such information on the quality of the correlation 
cannot be extracted from the q map. In order to 
investigate both the slope and the quality of the 
correlation, i.e. the correlation coefficient Tc, in SF 
regions and regions with no SF, we derive the clas- 
sical correlations in regimes of well-defined proper- 
ties. 

The radiation field is one of the properties which 
sets the dust heating. Since it is fed by UV photons 
from the SF regions as well as photons from old 
and low mass stars (that contribute significantly to 
the diffuse ISRF) , the radiation field could provide 
a proper defining parameter here. Taking advan- 
tage of the radiation field map (see Fig. 1, right), 
we derive the radio-FIR correlation separately for 
regions where the main heating source of the dust 
is either SF or the diffuse ISRF. Along the spiral 
arms and in the central disk, the radiation field has 
U > 1.3Uq, with Uq being the corresponding radia- 
tion field in the solar neighborhood, and it is smaller 
between the arms and in the outer disk (see Fig. 1). 
Hence, U — l.St/© is our criterior|3 to differentiate 



This is determined by considering the radiation field 
in the solar neighbourhood as a proxy to the diffuse 



regimes with ISRF-fed from SF-fed radiation fields. 
Table 7 lists the results of the 20 cm synchrotron- 
FIR correlation analysis for both monochromatic 
and bolometric FIR emissions in the two defined 
regimes. 

The slope b of the correlations is indeed shal- 
lower in the ISRF than that in the SF regime, in 
agreement with the q map. However, apart from 
the synchrotron correlation with the warm dust 
(70 /im) emission, the correlations in the ISRF 
regime are as good as those in the SF regime (com- 
pare Tc values). This is also the case for the bolo- 
metric FIR-synchrotron correlation. 

It is worth mentioning that, among the FIR 
bands, the 100 /zm emission has about the same 
slope as the bolometric FIR versus the synchrotron 
emission. This can be explained by the fact that the 
peak of the FIR SEP oc curs around 100 /xm in this 
galaxy (jPale et al.ll2012D . 

7.2. Correlations versus magnetic fields 

In the inter-arm regions of NGC 6946, the two mag- 
netic arms are well traced in the 20 cm synchrotron 
map (Fig. 3) and show lower q values. A lower q 
or a smaller FIR/synchrotron ratio indicates that 
the synchrotron emission is less influenced than the 
FIR emission by the absence of SF in the inter- 
arm regions. This can be explained by different de- 
pendencies of the FIR and synchrotron emission on 
SF, or that other mechanisms than SF regulate the 
synchrotron emission. The synchrotron emission is 
a function of the total magnetic field strength and 
CRE density. Thus, using the synchrotron emission 
as a tracer of SF implies that both the magnetic 
field and the CRE density are related to SF. This 
seems plausible since supernovae produce CREs 
and induce a turbulent ma gnetic field via strong 
shocks in SF regions (e.g. iRevnolds et"alll2012[ ). 
However, CREs diffuse and propagate away from 
their birth places into the magnetized ISM on large 
scales, where the total magnetic field is dominated 
by the uniform or ordered field. The origin of the 
ordered magnetic field can be li nked to the dy- 
namo effect on galactic scales (e.g. 'Beck et al."1 990l 
,1996.) and is not correlated with SfH ( c-g- iChvzvl 
I2OQ8I : lKrausell2009l: iFletcher et aI1l201l[ ). Therefore 
in NGC6946, magnetic arms seem to compensate 
for the lack of synchrotron emission in the inter- 
arms as shown by the anti-correlation between q 



ISRF optimized for enough number of points in be- 
tween the arms and outer disk of NGC 6946 at 70 nm 
(a 5(T level statistically significant correlation or t=ll 
for 70 /im-synchrotron correlation). 

^ The ordered magnetic field can also be produced by 
compressions and shear of the anisotropic or turbulent 
field in dense gas, known to be the origin of the strong 
ordered fiel d in the ce ntral 6kpc gas concentration of 
NGC 6946 l|Beckll200'7D 
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Table 7. Synchrotron-FIR/gas correlation in different regimes of radiation fields J/isrf and Usf- 



Y 




' C 








' c 






-i-70 

Iioo 
I160 
I250 


U.OOZIIU.U4: 

1.00±0.04 
0.90±0.02 
0.78±0.02 


n fii _i_n AC 

0.73±0.04 
0.79±0.02 
0.80±0.03 


1 1 n 
16.7 
27.4 
23.4 


245 
444 
396 


1 Qs-i-n 04 

1 .OOILU.U'4: 

1.32±0.04 
1.15±0.04 
1.07±0.03 


SO-t-O 09 

U.OUZIZU.UZi 

0.83±0.03 
0.85±0.02 
0.85±0.02 


9^^ A 

31.7 
35.7 
34.5 


449 
457 
450 


FIR 

Scas 


1.05±0.04 
0.68±0.02 


0.82±0.02 
0.59±0.03 


21.0 
18.6 


231 
648 


1.33±0.04 
1.00±0.03 


0.87±0.02 
0.86±0.02 


38.6 
36.1 


471 
481 


EsFR-Eoas 
Btot-Eoas 


1.45±0.06 
0.27±0.05 


0.58±0.04 
0.48±0.05 


14.69 
9.7 


320 
314 


1.34±0.05 
0.23±0.01 


0.64±0.04 
0.75±0.03 


17.8 
23.6 


452 
429 



Notes. Pearson correlation coefficients Vc and bisector slopes 6, log(Y)oc 61og(X), between the 20cm synchrotron 
(X) and various FIR bands, bolometric FIR, and total neutral gas (Y). The number of independent points, n, 
and the Student's t-test, t, are also indicated for each correlation. Also shown are the correlations between the 
gas surface density Ecas and the star formation surface density Espr, the total magnetic field Btot in logarithmic 
space. The nucleus was excluded before cross-correlating. 



and the linearly polarized intensity PI contours in 
Fig. 8 (middle panel). 

It is also instructive to investigate the role of 
the magnetic fields responsible for the synchrotron 
emission in places with high q (and high SF). 
Figures 6 and 8 show a good correspondence be- 
tween high q regions and Btur, Btot- The scatter 
plot in Fig. 9 (a) shows that q versus Btot and Btur 
obeys the following power-law relations: 



and 



goc (0.96 ±0.02) log Bt 



q(x (0.85 ±0.02) log Btur, 



(10) 
(11) 

obtained using the bisector method. Both fits 
have a dispersion of 0.07 dex. Hence, the FIR-to- 
synchrotron flux ratio linearly changes with the to- 
tal magnetic field strength. Studying similar linear- 
ity and also the flatter g-Btur correlation in other 
galaxies would be of high interest. 

7.3. Correlations versus star formation 

We now aim to characterize the relation between 
q and SF. We first derive the SF surface density 
(Ssfr) using the de-reddened Ha map. The dif- 
fuse emission is excluded by masking regions where 
the free- free 20 cm flux falls below a threshold value 
(0.05 mjy, that is th e minimum flux o f the de- 
tected Hii regions by iLacev et al.lll997l ). The nu- 
cleus was also subtracted as the optically thin con- 
dition does not hold there (rna > 1) and also be- 
cause of its anoma lous properties (for details see 
iMurphv et al.ll20TTl ). The de- reddened Ha luminos- 
ity is_com;ertedto SFR using the relation calibrated 
bv IMurphv eFall (|2011ll for NGC6946: 



SFR 



(Moyr-i; 



= 5.37 X 10" 



ergs ^ 



(12) 



where Lh^ is the H^ luminosity density. Ssfr (in 
units of Moyr^-'^ kpc^^) is then calculated. The 



power-law behaviour of q versus Ssfr is given by: 

goc (0.115 ±0.003) log Ssfr, (13) 

with a Pearson correlation coefficient of rc=0.8 and 
a dispersion of 0.05 dex (Fig. 9, b). Hence, the FIR- 
to-synchrotron fiux ratio changes with a much flat- 
ter power-law index with Ssfr than with Btot or 

Btur- 

7.4. Connection between star formation and 
magnetic fields 

As discussed above, star formation can regulate the 
magnetic field in galaxies. Therefore, understanding 
their relationship is important. In Fig. 9 (c), we 
show a good correlation between Ssfr and both 
Btot [re = 0.69 ±0.03) and Btur (re = 0.67 ±0.03). 
The bisector fit leads to a power-law proportionality 
for both Btot and Btur vs. Ssfr: 



log Btot a (0.14 ±0.01) log Ssfr, 



(14) 



log Btur (X (0.16 ±0.01) log Ssfr, (15) 

with dispersions of 0.08 and 0.06 dex, respectively. 
Eq. 14 reflects a process where the production of 
the total magnetic field scales with SF activity. This 
process probably has a common origin with the pro- 
duction of Btur as indicated by the similar slopes 
in Eqs. 14 and 15. Feedback mechanisms associ- 
ated with star formation, such as supernova and 
strong shocks, produce an increase in turbulence. 
These could amplify small-scale magnetic fields by 
a turbulent dynamo mechanism where kinetic en- 
ergy converts to m agnetic energy (e.g. iBeck et al.l 
FiggalGressel et al.ll200a . 

On the other hand, there is no correlation be- 
tween Bord (as traced by polarized intensity) and 
Ssfr for the whole range of parameters (tc = 
0.14 ± 0.02, see the crosses in Fig. 9d). However, 
a weak anti-correlation {tc = —0.25 ± 0.5) with a 
slope of 6 = —0.24 ±0.1 is found for high Ssfr 
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Fig. 9. The radio-IR correlation (q-parameter) versus total and turbulent magnetic field strength (a) 
and star formation surface density Ssfr (&)■ The magnetic field strength-SgpR relation is shown for the 
total, turbulent (c), and the ordered magnetic field strength measured by the degree of polarization (d). 
The lines in the upper panels are the bisector fits. The solid lines in c and d show the ordinary least 
square fit, and the dashed line the bisector fit. In panel d, the binned data points (diamonds) are fitted 
for SsFR ^ 0.01 Mq yr^^ kpc^^ (solid line). In panels a and c, the total magnetic field is shifted up by 
0.3 units for clarity. The nucleus is not included in these analysis. 



(> 0.01 M0 yr^^ kpc^^), after binning the data 
points with a width of 0.5 (see the red diamonds in 
Fig. 9d). An anti-correlation for large Ssfr is also 
visible when comparing the PI and the SF contours 
in Fig. 8, where no polarized emission is associated 
with the optical spiral arms. Using the wavelet anal- 
ysis, an anti-correlation bet ween PI and Haemis- 
sion was already found by iFrick et "all ()2001[ ) on 
scales equivalent to the width of the spiral arms. 
All taken together, this suggests that the produc- 
tion of Bord could be suppressed along the spiral 
arms and in SF regions. An efficient dynamo action 
in the inter-arm regions could produce a regular 
field that is anti-p roportional to Ssfr,. as show n 
for NGC6946 (see lBeca[2007HRohde et al.lll"999h . 

Comparison with a similar study in another spi- 
ral galaxy, NGC 4254, is instructive. This galaxy is 
a member of the Virgo galaxy cluster that is ex- 
periencing a gravit ational encounter at the clus- 
ter's periphery (e.g. lChvzv et al.ir2007D or ram pres- 



sure due to the ga laxy motion throug h the in- 
tracluster medium (Murphy et al.l l2009l) or both 
of them (iVoUmer et al. 2005). Both effects, tidal 
forces and ram pressures can compress or shear 
the magnetic fields. Interestingly, NGC 4254 also 
shows two different trends of Bord v s. SsfRi wit h 
no general correlation between them (jChvzvll2008[ ). 
The division of the two-way behavior, however, oc- 
curs at a larger Ssfr (— 0.02 yr"^ kpc~^) than 
in NGC 6946. Moreover, in NGC 4254, the slope 
of Btur-S^SFR relation (0.26±0.01) is larger than 
that of Btot-SsFR relation (0.18±0.01), unlike in 
NGC 6946. The slope of Btur-SsFR relation is also 
larger in NGC 4254 than in NGC 6946 (0.16±0.01, 
see Eq. 15). Hence, although the two galaxies re- 
semble each other in their general magnetic field- 
SsFR behavior, they differ in the details of these 
connections. These differences could be related to 
different contributions of Bord and Btm- to Btot, 
as a similar scaling relation holds between Btot 
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and EsFR in these galaxies. Hence, investingat- 
ing Bord-to-Btur ratio (Bord/Btur) could be instruc- 
tive. In NGC4254, Bord/Btur shows a general de- 
creasing trend with increasing S spp for the w hole 
range of values (see Fig. 7b of IChvzvl I2008D . In 
NGC 6946, such a decreasing trend is indicated only 
for EsFR > O.OlMQyr-i kpc'^ (Fig. 10). For smah 
5^SFR values, Bord/Btur IS larger in NGC 4254 than 
in NGC 6946. 

In the Virgo Cluster galaxy NGC 4254, the 
strongest ordered field is found in the outer arms, 
dominated by anisotropic turbulent field and is 
probably a product of shearing /stretching forces 
cause d by weak gravitatio nal interaction . (jChyzvl 
|2008[1 and/or ram pressure (jMurphv et al.ll2009D in 
the cluster environment. Such forces can transform 
Btur into Bord in the outer arms of NGC 4254 where 
SsFR is low (particularly in the southern arm) . This 
could cause the steeper log(Btur)-log(SsFR) than 
log(Btot)-log(EsFR) relation in NGC 4254. 

In spite of the magnetic arms, NGC6946 yields 
smaller values of Bord/Btur at low SFR compared 
to NGC4254 because it does not experience shear- 
ing due to tidal forces and/or ram pressures in- 
serted from a cluster environment. The total field 
has similar properties in both galaxies (same slope), 
because the shearing does not change the total 
field strength. It only transforms turbulent into 
ordered field. This could explain the flatter slope 
of log(Btur)-log(i;sFR) in NGC6946 than that in 
NGC4254. 

Therefore, we might expect that, in normal 
galaxies like NGC6946, the slope of log(B) vs. 
log(I]sFR) is similar for both Btur and Btot, while it 
is steeper for Btur than for Btot in galaxies interact- 
ing with the cluster environment (like NGC4254). 



7.5. Correlations with neutral gas 

Some of the radio-FIR correlation models are based 
on a coupling between the magnetic field B and 
gas d ensity p (H elou fc B icav"1993': 'Niklas & BecW 
Il997t iHocrncs et al.l 11998: Murgia et al. 2005). 
Physically, the B-p coupling could be caused by 
the amplification of the Magneto Hydrodynamic 
(MHD) turbulence until en ergy equipartition is 
reached (jCroves et al.l[2003h . This model has the 
advantage of 1) producing the radio-FIR correla- 
tion where the FIR is dominated by dust heated 
by older stellar po pulations or the ISRF (e.g. 
iHoernes et al.l [19981 ) and 2) naturally producing a 
local radio-FIR correlation, breaking down only on 
scales of the diffusion length of t he CREs Idu (as 
discussed in iNiklas fc BecMll997D . 

Table 7 shows that the total molecular and 
atomic gas surface density Scas in Mq pc~^ is cor- 
related with the synchrotron emission in both the 
SF and the ISRF regimes. The correlation, in the 
ISRF regime, is however not as tight as in the SF 



regime (see Tc values in Table 7) because the inter- 
arm regions (heated by the ISRF) are filled with 
low-density atomic HI gas which appears to have a 
weaker correlation with synchrotron emission com- 
pared to dens e mole cular gas (Beck 2007,). In M51, 
iTilanus et al.l ()1988D found a much weaker correla- 
tion between the diffuse HI gas and the synchrotron 
emission. 

A tight correlation between Scas and the total 
magnetic field holds for the entire galaxy (Fig. 11). 
The total magnetic field Btot vs. Ecas follows a 
power-law relation, 

(16) 



log Btot (X (0.23 ±0.01) log Ecas. 



The dispersion of the fit is about 0.04. In the ISRF 
regime, the Btot~EGas correlation is weaker than in 
the SF regime, although the power-law indices are 
si milar (Table 7). 

INiklas fc Bec3 (|1997D suggested that the tight 
radio-FIR correlation can be reached by consid- 
ering the coupling between the FIR emission and 
the gas density via the Kennicutt - Schmidt (KS) 
relation be tween Esfr and Scas (jSchmidtl Il959t 
lKennicuttl [T998). Interestingly, in NGC 6946, not 
only does the KS relation hold, but also it is very 
similar in the ISRF and SF regimes (Table 7 and 
Fig. 12). Hence, it seems that in the ISRF regime, 
the lower gas density and the weaker radiation 
field conspire to hold the KS relation but with a 
weaker EsFR~5]Gas correlation. In the SF regime, 
the EsFR- Egas power-law in dex of ~ 1.3 is in agree- 
ment with lKennicuttI (|1998D . The KS index in both 
regimes agrees with the gl obal value o f 1.46i b0.29 
for NGC 6946 presented in lBigiel et aP (120081) . 

An equipartition between the CRE and mag- 
netic fiel d energy densi ti es is needed in the B-p 
model of INiklas fc Beckl (|1997D . This implies that 
the energy density of CREs (and the synchrotron 
luminosity) is determined by the field strength 
alone, and neither the CRE production rate nor the 
CRE escape probability affects it. On small spatial 
scales, however, e.g. in supernova remnants, the en- 
ergy density of the particles may exceed strongly 
the energy density of the magnetic field. A static 
pressure equilibrium could then be achieved when 
the components of the interstellar medium have 
been relaxed on larger scales (e.g. determined by 
a scale height of gas clouds of a few 100 pc). Hence, 
the local synchrotron-FIR correlation should break 
down on small spatial scales where equipartition 
is not valid anymore. Looking for such a break 
in the synchrotron-FIR correlation in NCG6946, 
we refer to our scale-by-scale correlation study. 
Figure 7 shows that the smallest scale on which the 
synchrotron-70 /zm correlation holds, i.e. r^, ~ 0.5, 
is a ~ 1.7 kpc (on scales smaller, r„ < 0.5 and hence 
it is below the threshold value for the acceptabil- 
ity of the correlation) . This scale can be translated 
as ^dif or the scale of the static p ressure equilib- 
rium, according to the B-p model (|Niklas fc Beckl 
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Fig. 12. Correlations between the Ssfr [M© yr~^ kpc~^] and the total gas surface density Scas [M© pc~^] 
for the ISRF (left) and SF regimes (right) as defined by /7 > \.?>Uq. Also shown are the bisector fits. 
The nucleus is not included in this analysis. 



Il997f ). However, this is not certain, since the colder 
dust traced by the longer FIR wavelengths do 
not show similar break r^, < 0.5 when correlat- 
ing with the synchrotron emission. Nevertheless, 
the synchrotron-70 /im correlation has been used 
in previous studies to determine the CRE diffu- 
sion length. For instance, using the Spitzer MIPS 
70 ^m data, iMurphv et al.1 (|2008f ) estimated ^dif in 
a sample of nearby galaxies by applying an image- 
smearing method. Using the same FIR wavelength 
(70 /zm) , our derived ^dif agre e s with the best- fitted 
disk value of iMurphv eFaD (|2008[ ) in NGC6946, 
^dif = 1-6 ± 0.1 kpc, although the two methods are 
different. 



7.6. The origin of the radio-FIR correlation in 
NGC6946 

The proposed models to explain the radio-FIR cor- 
relation in galaxies include: 

a) The calorimeter model of IVolkl (|1989D . This 
model assumes that CREs are completely trapped 
in their host galaxies and that the galaxies are op- 
tically thick for the dust-heating stellar UV pho- 
tons. In the case of complete CRE trapping, a syn- 
chrotron spectral index of a„ > 1.0 is expected, 
i.e. synchrotron and inverse Compton losses are the 
dominant cooling mechanism of CREs. 

b) Th e o ptically-thin in odel of iHelou fc Bicavl 
(|l993h and lNiklas fc Beckl (I1997D . The dust-heating 
UV photons and the cosmic rays are assumed 
to have a common origin in m assive star forma- 
tion similar to the model of IVolkl ([l989), but 
in this case both photons and cosmic rays can 
escape the galaxy. Although different distribu- 
tion for the CREs e nergy density is app lied b; 



IHelou fc Bicavl (I1993D and 'Niklas & BecW ^199 
(also see lHoernes et a l. 1998; Murgia et al. 20051 
the basic condition of B-p coupling is a common as- 
sumption in these models, as di scussed in Se c t. 7.5 . 
c) The calorimeter model of iLacki et al.l (|2010( ) 



in starbursts, in which cosmic ray protons lose 
all of their energy and produce secondary elec- 
trons and positrons whose synchrotron emission 
keeps the radi o-FIR correlation linear. This was 
motivated by 'Thompso n et al.l ([2006) who sug- 
gested that Bremsstrahlung and ionization losses 
are more important in starburs t galaxies leading t o 
a flat an (~0.75). According to lLacki etaP (|2010[ ). 
the CRE density is not proportional to the SFR 
in starbursts . in contras t to both the calorime- 
ter mode l of Volk |l989l) and the 'optically-thin' 
model s of lHelou fc Bicavl (|1993D and lNiklas fc Beckl 
(|1997f ). who explained the observed non-linear 
synchrotron-FIR corre lation in the Shapley-Ames 
galaxies (presented in lNiklasl[l997a^ by assuming 
a SFR related synchrotron emission and B-p cou- 
pling. 

In NGC6946, the complete CRE trapping as- 
sumed in the model (a) does not apply, since de- 
pending on the location, the synchrotron spectral 
index changes between a flatter a„ ~0.6 (in the 
star-forming regions) and a steeper ~ 1 (in re- 
gions of strong ordered magnetic field. Sect. 4). This 
indicates the presence of various loss mechanisms 
of the CREs energy leading to a flatter nonther- 
mal spectrum (mean an ~0.8) than that produced 
by the efficient synchrotron and inverse Compton 
losses (q;„= 1.0-1.2) for the entire galaxy. 

Although NGC 6946 is classified as a normal 
star-forming galaxy, it would be instructive to 
see its position in respect to the model (c). The 
synchrotron-FIR correlation is non-linear with the 
warmer dust globally (Table 5), or the dust heated 
by the SF (Table 7). On the contrary, a linearity is 
found for the colder dust globally (Table 5), or the 
dust heated by the ISRF (Table 7). The linearity of 
the correlation cannot be explained by secondary 
electrons fc positrons suggested for dense gas con- 
dition in starbursts, due to the low gas density of 
the ISRF regi me in NGC 6946 (in agreement with 
lMurphvll200l . 
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As shown in Sect. 7.5, we already obtained indi- 
cations for a B-jO coupling in NGC 6946 which is in 
favor of the model (b). How this model could also 
reproduce the observed slope of the synchrotron- 
FIR correlation in the SF and the ISRF regimes 
are detailed as follows: 

In the SF regime, the KS relation gives 
SsFR oc Sg^,, with the KS index of 1.3±0.05 
(Table 7). Our data show that the FIR luminosity 
is proportional to Ssfr with a slope of 0.95±0.05. 
Hence, 

FIR oc 

We note that the same proportionality is derived by 
directly correlating between the FIR luminosity and 
5^SFR, as expected. From Eq. (6), the equipartition 
magnetic field is given by: 



with an 
that 



B, , CX /l/(3+a„) 
■^tot ^ -'syn J 

0.8±0.1 (Sect. 4). In Sect. 8.2, we found 



Btnt oc E 



0.23±0.01 

Gas 



Follow ing iNiklas fc Bec5 (|1997D and lHoernes et al.l 
(|1998() . we assume that the scale height of the neu- 
tral gas is equal to that of the dust and is constant. 
Thus, Scas can be replaced by the gas volume den- 
sity in the above proportionalities resulting in 



FIR oc I^^t'^'-'\ 



(17) 



This is in excellent agreement with the observed 
synchrotron-FIR correlation in the SF regime with 
the slope of b^^ = 1.33 ±0.04 (Table 7). 

A similar calculation for the ISRF regime leads 

to 

(18) 



FIR oc 



The slope is higher than the observed slope of 
b^ ^ = 1.05 ± 0.04 by about 2a errors (28%). 

IXu et al.l (|1994l ) found a linear correlation be- 
tween the 20 cm synchrotron and cold dust emis- 
sion in a sample of late- type spirals. They explained 
this correlation by considering intermediate mass 
stars (5-20 Mq) as a heating source of cold dust 
and synchrotron emission, since these stars are su- 
pernova progenitor as well. It is also likely that, 
in NGC 6946 which is a late-type spiral with a 
good 20 cm synchrotron-cold dust correlation (e.g. 
see Table 5), the intermediate mass stars provide 
the non-ionizing UV photons to heat the cold dust 
(which emit at longer FIR wavelengths). As such, 
these UV photons provide the bulk of the ISRF in 
NGC 6946. 

In the previous sections, we showed that the lo- 
cal radio-FIR correlation varies as a function of 
not only star formation rate but also dust heat- 
ing sources, magnetic fields, and gas density; i.e., 
the ISM properties. The SF and ISM are contin- 
uously influencing each other: Stars form within 
the dense and cold regions of the ISM, molecular 



clouds, and replenish the ISM with matter and en- 
ergy. Hence, the local radio-FIR correlation is a 
probe for the SF-ISM interplay which is partly ex- 
plained through scaling relations such as the KS 
relation between Ssfr and gas surface density and 
the relation between Ssfr and the magnetic field 
strengths (Sect. 7.4). 

The global radio-FIR correlation is known to 
be a tracer of SF. Among the SF-based theories 
to explain the global radio-FIR correlation, our lo- 
cal studies match with those considering a coupling 
between the magnetic field and gas density, as the 
radio-FIR correlation also holds in regions with no 
massive SF, e.g. in the inter-arm regions and the 
outer disk (Sect. 7.1). This shows that a balance be- 
tween the gas and magnetic field/CRE pressure is 
an unavoidable condition for the correlation. Apart 
from this, the global radio-FIR correlation as a 
tracer of SF still applies since the integrated radio 
and FIR fluxes are weighted towards more luminous 
regions of galaxies, i.e. the SF regions. 

7.7. ISM and propagation of CREs 

Propagating through the ISM, CREs can ex- 
perience various energy losses via ionization, 
bremsstrahlung, adiabatic, synchrotron, and 
inverse-Compton losses that change the power-law 
index of the energy distribution of these particles 
or equivalently the nonthermal spectral index, 
an. The maps of the synchrotron spectral index 
(Sect. 4), magnetic fields (Sect. 5), and radiation 
field (Sect. 3), provide direct information on the 
main cooling mechanisms of CREs, as well as 
the cooling timescale and diffusion scalelength of 
CREs in NGC 6946. As shown in Fig. 4, there is 
a difference in the spectral indices of Aa^ ~0.5 
between the star-forming regions and the magnetic 
arms/inter-arm regions. This is expected if the 
main mechanisms of energy losses for the electrons 
are synchro tron em i ssion and inverse Compton 
scattering ()Longaiij 119941) . The CRE cooling 
timescale icooi (in units of yr) associated with the 
two dominant pro cesses can be deriv ed from the 
following formula (|Murphv et al.ll2008[ ): 



-0.5 



cool 



5.7 X 10 



7 (J^X 

Vghz/ 

Ub + C/rad 



>G^ HO-12 



(19) 



■ erg cm 

where ~ BE^ is the critical frequency at which 
a CRE emits most of its energy E. Ub = B^/(87r) 
is the magnetic field energy density, and C/rad is the 
radiation energy density. 

Using a simple random walk equation, CREs 
will diffuse over a distance ^cooi = (-Db tcooi)°'^ be- 
fore losing all of their energy to synchrotron and 
inverse Compton losses, with the energy-dependent 
diffusion coefficient De- Assuming that the dif- 
fusion length scale of 1.7 kpc, obtained from the 
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synchrotron-FIR correlation (Sect. 7.5), is equiva- 
lent to the cooling scale length 4ooii De can be es- 
timated independently. iMathis et al.l ()1983[ ) deter- 
mined the radiation field energy density in starlight 
of J7rad = 8.64 X 10^^^ J7 erg cm~^ . However, 
f/rad should also include the dust emission en- 
ergy density, Ufyr. = 5 x 10~^^ J7 ergcm"-^ as 
well as the cosmic microwave background radia- 
tion ene rgy density, Ucmb — 4.17 x 10"^'^ ergcm^"^ 
(|Drainel[20in . Therefore, we use a more general 
formof{7rad = (4.17 + 13.64 [/) X 10~"ergcm-3. 
For the entire galaxy, U ~2 [/© and B ~ 18 fxG. 
Thus, Ub ^ 1.3 X 10-"ergcm-3 and t/^ad ^ 
3.1 X 10^^^ ergcm"'^. Substituting these values 
in Eq. (19) results in tcooi ^ 1-3 x 10^ yr. The 
CRE diffusion coefficient is then De — 6.8 x 
lO^^cm^s"^. Assuming that De changes w ith en- 
ergy as De = Do {E/GeV)'^-^ (e.g. Ginzburg " et al.l 
llQSOh . we derive the normalization factor Dq = 
4.6 X lO^^cm^s-i for the 2.2 GeV CREs. 

As discussed in Sect. 6.2, the drop in the 
synchrotron-FIR correlation on the larger scales of 
about 2.5 kpc is related to a lack of FIR emission 
from the magnetic arms (in the inter- arm regions), 
where diffused CREs experience synchrotron loss. 
Hence, we assume that the cooling length scale is 
?cooi — 2.5 kpc in the magnetic arms, where the ra- 
diation field is ?7 ~ Uq and the magnetic field 
B~ 15 /iG. These lead to a cooling time scale of 
icooi = 1.7 X 10^ yr and a diffusion coefficient of 
Do = 7.0 X 10^8 cm^ s'^ for the 2.4 GeV CREs. 

The above estimates of the CRE diffu- 
sion coefficient agree with both observationa l 
([Strong fc Moskale nkol 119981; 'Dahlem |t all I1995D 
and theoretical (e.g. Roediger et al. 200^ estimates 
of ~ (1 - 10) X lO^^cm^s-i and is about 10 
times larger than the diffusion coefficient estimated 
on small scales in the turbulent medium near 30Dor 
([Murphv et al .1120121 ). These demonstrate the effect 
of the ISM and CRE cooling on the synchrotron- 
FIR correlation on kpc-scales. On these (i.e., > 
Ikpc) scales, the CRE population is likely to be 
dominated by old/diffuse particles accelerated in 
past generations of star formation. These parti- 
cles may even be re-accelerated by passing shocks 
in the ISM, and thus will have little memory of 
their original birth sites in star-forming regions. 
Consequently, the propagation of these particles is 
very sensitive to the ISM conditions. 

On sub-kpc scales, however, the CRE popula- 
tion is more likely to be dominated by younger 
CREs, which are still close to their production sites 
in star-forming complexes. The energy distribution 
of these particles is more influenced by the age and 
activity of their acceleration sources rather than 
the quasi-state ISM condition. This is in agreement 
with Murphy et al. (2008) who showed that the 
current diffusion length of CREs from star-forming 
structures is largely set by the age of the star- 



formation activity rather than the cooling mech- 
anisms in the general ISM. 

8. Summary 

Highly resolved and sensitive Herschel images of 
NGC6946 at 70, 100, 160, and 250 ^im enable us to 
study the radio-FIR correlation, its variations and 
dependencies on star formation and ISM properties 
across the galactic disk. Our study includes differ- 
ent thermal/nonthcrmal separation methods. The 
radio-FIR correlation is calculated using the classi- 
cal pixel-by-pixel correlation, wavelet scale-by-scale 
correlation, and the q-method. The most important 
findings of this study are summarized as follows: 

- The slope of the radio-FIR correlation across 
the galaxy varies as a function of both star for- 
mation rate density and magnetic field strength. 
The total and turbulent magnetic field strengths 
are correlated with SgpR with a power-law in- 
dex of 0.14 and 0.16, respectively (Eqs.l4, 15). 
This indicates efficient production of turbulent 
magnetic fields with increasing turbulence in ac- 
tively starforming regions, in general agreement 
with Chvzv (2008). 

- In regions where the main heating source of 
dust is the general ISRF, the synchrotron emis- 
sion correlates better with the cold dust than 
with the warm dust. However, there is no differ- 
ence between the quality of the correlations for 
colder /warmer dust in regions of a strong radi- 
ation field powered by massive stars. This is ex- 
pected if warmer dust is mainly heated by SF re- 
gions (where synchrotron emission is produced 
by young CREs/turbulent magnetic field), and 
colder dust by a diffuse ISRF across the disk 
(where synchrotron emission produced by old, 
diffused CREs/large-scale magnetic field). 

- The synchrotron-FIR correlation in strong ra- 
diation fields can be well explained by the 
optically-thin models where massive stars are 
the common source of the radio and FIR 
emission. The intermediate-mass stars seem 
to be a more appropriate origin for the ob- 
served synchrotron-FIR correlation in the ISRF 
regime. 

- The synchrotron spectral index map indicates a 
change in the cooling of CREs when they propa- 
gate from their place of birth in star-forming re- 
gions across the disk of NGC 6946. Young CREs 
emitting synchrotron emission with a flat spec- 
trum, an — 0.6 ± 0.1, are found in star-forming 
regions. Diffused and older CREs (with lower 
energies) emit synchrotron emission with a steep 
spectrum, q;„ = 1.0 ± 0.1, along the so-called 
'magnetic arms' (indicating strong synchrotron 
losses) in the inter-arm regions. The mean syn- 
chrotron spectral index is a„ — 0.8 ± 0.1 across 
the disk of NGC 6946. 
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- The cooling scale length of CREs determined us- 
ing the multi-scale analysis of the synchrotron- 
FIR correlation provides an independent mea- 
sure of the CRE diffusion coefficient. Our deter- 
mined value of Do = 4.6 x 10^^ cm^ for the 
2.2 GcV CREs agrees with the observed values 
in the Milky Way. This agreement suggests that, 
reversing our argument and assuming Milky 
Way values for De, the cooling scale length 
of CREs due to the synchrotron and inverse- 
Compton energy losses appear to be consistent 
with scales on which the radio-FIR correlation 
is weak on kpc scales. This indicates that the 
interstellar magnetic fields can affect the prop- 
agation of the old/diffuse CREs on large scales. 
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radio frequencies and assumes a priori knowledge of . An 
incomplete frequency coverage of the radio images makes it 
difficult to find variations in «„, technically, and hence in 
most cases the thermal/nonthermal components are decom- 
posed assuming a constant an- Although this assumption 
does not allow for the study of the change in the CREs en- 
ergy when they propagate away from SF regions, the stan- 
dard separation method is often used as it can produce ther- 
mal and nonthermal maps using radio continuum data at 
only two frequencies (at high and low radio band) and more 
straightforwardly than using, e.g., the TRT method. Here, 
we derive the thermal/nonthermal distributions from the 
standard method for NGC6946. Comparing the results of 
the standard and the TRT methods could be instructive for 
similar studies in future. 

For a total spectral index, a, obtained from the observed 
flux densities at frequencies ui and U2 and constant value of 
the nonthermal spectral index an, the thermal fraction at 
frequency ui is given by 



Ul VI Ul Ul 

ein et all Il984f) . Then the thermal flux density at fre- 



the nonthermal flux density S^^ 



S-^ X and 



Using the 



data at 3.5 and 20 cm for the above formula, and assum- 
ing the freque ntly adopt ed index a„ — 1 (e.g. iKlein et al.l 
1984; Berkhui jsen et al.ll2003i ). the correspon ding free-fre e 
and synchrotron maps are derived (see also iBeckl l2007f l . 
Figure A.l shows the resulting maps at 3.5 cm divided by 
those obtained in Sect. 3 (the standard-to-TRT ratio maps). 
The standard-to-TRT free-free ratio fluctuates around one 
in the spiral arms and the central disk indicated by the TIR 
contours. In the nucleus and its surroundings, however, the 
ratio exceeds 2. A striking excess of free-free emission from 
the standard method occurs between the arms in the west 
and the south. We note that, in these regions, there is no 
significant emission of Ha as a tracer of the ionized gas, nor 
24 or 70 fim as tracers of warm dust, and even from the 
TIR as a tracer for SF. Hence, the strong diffuse free-free 
emission in those regions is not real. 

The standard-to-TRT synchrotron ratio map resembles 
the an map (Fig. A.l). It is close to unity in regions where 
the assumed an is close to that obtained from the TRT 
model. The synchrotron emission is underestimated in re- 
gions with a„ flatter than the assumed a„ = 1. The lack of 
synchrotron emission is particularly seen in regions with an 
excess of free- free emission in the west /south of the galaxy 
discussed above. The synchrotron emission from the SF re- 
gions and the nucleus is also underestimated as < 1 for 
those regions (Sect. 4). 



Appendix A: Thermal/nonthermal 

separation using the standard method 

iTabatabaei et aP ll2007bfl presented a detailed comparison of 
the TRT thermal/nonthermal separation method to a more 
conventional method, the so called 'standard method', for 
M33. The standard method requires accurate absolute mea- 
surements of the brightness temperature at high and low 
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Fig. 7. The scale-by-scale correlation between the 
FIR emission and the observed 20cm RC (top), 
the free- free (middle), and the synchrotron (bot- 
tom) emissions in NGC6946. The horizontal line 
shows — 0.5, the threshold for detecting a cor- 
relation. Larger r^ia) corresponds to more highly 
correlated structures. These figures illustrate that, 
on scales smaller than 4 kpc, the FIR-synchrotron 
correlation is weaker than the FIR-free-free corre- 
lation. 
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Fig. 8. The synchrotron-FIR correlation (q) over- 
laid with contours of the star formation traced by 
free-free emission with contour levels of 20, 50, 
100, 200, 400/^Jy/beam (top), the ordered mag- 
netic field traced by the linearly polarized intensity 
PI with contour levels of 85, 120, 160/j.Jy/beam 
(middle) and the turbulent magnetic field (/syn — 
PJ/0.75) with contour levels of 0.3, 0.6, 1, 2, 25, 
32mJy/beam (bottom) in NGC6946. 
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Fig. 10. The ratio of the Ordered-to-turbulent 
magnetic field strengths Boid/Btur vs. star for- 
mation surface density Esfr in NGC6946. The 
dashed hne shows the corresponding fitted hne in 
NGC4254 shifted by 0.2 a long the x-a xis {y = 
-0.32 a; -0.7, see Fig. 7b in lChvzvl[2008h . 



1.5 



3 1.3 





















^i**^* r = 0.80+-0.02 







0.5 1.0 1.5 2.0 
LOC Eg^g [Ma pc- 



2.5 



3,0 




Fig. A.l. Top: the ratio map of the standard-to-TRT 
based free-free intensity overlaid with contours of the TIR 
emission. The contour levels are 9, 15, 20, 130, 200 Lq kpc~^. 
Bottom: the same for the synchrotron component. The ratio 
values are indicated by the bars at the right of each image. 



Fig. 11. Correlation between the magnetic field 
and the total gas surface density Scas- Also shown 
is the bisector fit. The nucleus is not included in 
this analysis. 
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